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Abstract 

We  develop  a  class  of  multivariate  life  distributions  based  on  the  notion  that 
an  unknown  common  environment  induces  dependencies.  The  environment  is  assumed 
to  be  dynamic  and  is  described  by  a  gamma  process.  The  multivariate  exponential 
distribution  of  Marshall  and  Olkin,  motivated  via  an  environment  comprising  of 
nonfatal  Poisson  shock  processes  wherein  the  probability  of  failure  associated  with 
each  shock  is  the  same,  turns  out  to  be  a  special  case  of  our  development.  The 
gamma  process  environment  reduces  to  a  consideration  of  a  Poisson  shock  model 
with  the  shocks  having  different  probabilities  of  failure. 
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1 .  Introduction 


Consider  a  two  component  system  which  operates  in  an  environment  which 
may  or  may  not  be  identical  to  the  test  bench  environment.  The  component 
lifelengths,  when  assessed  under  the  test  bench  environment,  are  assumed  to  have 
failure  rates  \^(u),  i  =  l,  2,  for  u  ^  0.  The  operating  environment  is  supposed  to  be 
made  up  of  several  stresses  whose  intensities  and  presence  change  over  time.  It  is 
assumed  that  the  nett  effect  of  the  operating  environment  is  to  modulate  X-(u)  to 
X-(u)  77(u),  i  =  l,  2.  The  \^(u)’s  may  be  known  or  unknown;  however,  ?7(u)  is  assumed 
to  be  unknown  for  all  u  ^  0.  If  at  any  time  u,  the  operating  environment  is  judged 
harsher  (gentler)  than  the  test  bench  environment,  then  77(u)  >(<)!;  Viu)  =  1 
would  correspond  to  the  case  in  which  the  operating  environment  and  the  test  bench 
environments  are  identical.  Given  \j^(u),  X->(u)  and  n(u),  for  all  u  ^  0,  we  shall  judge 
the  component  lifelengths  independent.  When  ’7(u),  u  ^  0,  is  unknown,  we  shall 
describe  our  uncertainty  about  it  by  a  suitable  stochastic  process  (77(u);  u  ^  0}. 
Under  the  above  circumstances,  there  is  an  induced  dependence  between  the 
lifetimes  of  the  two  components,  and  the  aim  of  this  paper  is  to  describe  such 
lifelengths  for  a  particular  choice  of  {77(u);  u  ^  0},  namely,  that  {Jq  ?7(u)  du;  t  ^  0} 
is  a  gamma  process.  In  Singpurwalla  and  Youngren  (1989b)  we  consider  the  case 
wherein  {7?(u);  u  ^  0}  is  described  by  a  ”shot-noise  process”  [cf.  Cox  and  Isham 
(1980),  p.  34).  The  case  Xj(u)  =  X-,  i=l,  2,  with  r](u)  =  V,  and  uncertainty  about  V 
described  by  a  gamma  distribution,  has  been  considered  by  Lindley  and  Singpurwalla 
(1986),  Currit  and  Singpurwalla  (1989)  and  the  references  therein.  A  motivation  for 
describing  {Jq  ^(u)du;  t  ^  0}  by  a  gamma  process  is  in  Singpurwalla  and  Youngren 
(1989a),  and  in  Kalbfleisch  and  Prentice  (1980),  p.  203.  By  way  of  a  brief  review, 
and  with  the  aim  of  introducing  some  notation,  we  have 


Definition  1.1.  Let  a(t)  be  a  nondecreasing  ieft-continuous  real  valued  function  on 
[0,  ^),  with  0.(0)  =  0,  and  let  6  €  (0,  ^).  A  stochastic  process  (x(t);  t  ^  0)  is  said 
to  be  a  gamma  process  with  parameters  a(t)  and  B,  denoted  "XCt)  -  Q  (a(t),  5)  if 

i)  X(0)  =  0 

ii)  X(t)  has  independent  increments,  and 

jjj)  (X(t)  -  X(s))  has  a  gamma  distribution  with  shape  parameter 
(a(t)  —  a(s)),  and  scale  parameter  1/6,  for  0  s  ^  t. 


Definition  1.2.  [Dykstra  and  Laud  (1981)].  Let  Six),  t  ^  0,  be  a  positive  right- 
continuous  real  valued  function  and  let  X(t)  G  Q  (a(t),  1).  Then  the  process 

ft 

{z(t);  t  >  0),  with  z(t)  3(s)dx(s)  is  called  an  extended  gamma  process  with 

parameters  o(t)  and  .i3(t),  and  is  denoted  Z(t)  —  y£(o(t),  3i,tj. 

A  useful  property  of  the  extended  gamma  process  is  that  if  Z(t)  G  (o'A), 
5(t)),  then  G2(t)(w),  the  Laplace  Stieltjes  transform  of  the  distribution  of  Z(t),  is 
given  as 


(1.1) 


rt 

-  ( 

JO 


Cn  (1  —  w  ,3(u))  d  a(u) 


see  Dykstra  and  Laud  (1981)  or  Cinlar  (1980)  [who  introduces  another  generalization 
of  the  gamma  process  and  describes  an  imaginative  modeling  of  the  deformation  laws 
of  materials]* 


Let  Y(t) 


n(u)du,  and  suppose  that  Y(t)  G  g(a(t),  1/b),  where  a.(t)  is 


def 


continuously  differentiable  with  ^  a(t)  =  a(t).  Then,  it  follows  that  (t)  - 

ft 

1  X-(u)  d  Y(u),  the  cumulative  hazard  rate  of  component  i,  i==l,  is  such  that 
A^(t)  G  Qg  (a(t),  and  this  will  be  the  underlying  hypothesis  upon  which  our 


results  will  be  based* 


2.  The  Bivariate  Survival  Function  and  Its  Marginals 


Let  and  T-,  be  the  lifelengths  of  the  two  components,  and  suppose  that 
0  <  T  <  <  oo  .  Then,  under  the  hypothesis  of  the  previous  section  and  from 

the  independent  increments  property  of  the  gamma  process,  it  can  be  seen  that  the 
bivariate  survival  function  of  the  two  components 


def 


(2.1)  P(T^  >  Tp  T2  >  T2iX^(u),  '’'2)  ~ 


exp 


f  ^  £n  [1  +  j  a(u)du  •  exp|—  \  £n[l  +  -^)a(u)du 

Jq  ^  ^  •'^1 


also,  the  marginal  survival  functions  PCT-  >  T|Xj(u),  a(u),  b)  Fj(/) 

exp|^ —  f  £n(l  +  — )  a(u)duj,  i=l,  2. 

jQ 


2.1  The  Bivariate  Exponential  and  )X  eibull  as  Special  Cases 

In  what  follows,  we  show  that  one  of  the  most  widely  discussed  multivariate 
distribution  in  reliability  theory  -  which  has  been  motivated  via  considerations 
which,  at  least  at  the  surface,  appear  to  be  different  from  ours  -  is  a  special  case 
of  (2.1).  To  see  this,  suppose  that  X-(u)  =  X-,  i  — 1,  2,  and  all  u  ^  0;  then  Aj(  / ) 
X-Y(r)  £  q[cc(t),  ^],  from  which  it  can  be  verified  that 


(2.2) 

and  that 
(2.3) 


F(Tp  T 


:'  =  ( 


b+(X  j^+X2) 


ja(Ti)  j 


^  .a(T.^)-a(Tp 

b+X^J 


(xir) 


3 


Note  that  the  failure  rate  of  T^, 


(2.4)  r-Cr)  —  aCr)  £n 

is  a  constant  times  the  derivative  of  a(t),  the  shape  parameter  of  the  gamma 
process.  Thus,  we  may  choose  different  functional  forms  for  a(t),  each  suggested 
by  a  physical  scenario  of  the  environment,  and  obtain  different  marginal 
distributions  for  the  component  lifelengths.  One  such  choice  -  perhaps  a  natural 
one  -  is  that  o.(.r)  =  ar,  for  some  constant  a  >  0.  For  this  choice  of  aCr),  the 
marginal  distributions  are  exponential,  and  the  bivariate  survival  function  (2.2)  is 


^  b+X-i  .|0cT ^ 

(  ^+^1  ' 

jCXT^  ^  b(b+X|+X^)  ^^^^2 

^b+X 

tb+Xi+X^' 
1  ^ 

1  Hb+X^)(b+XpJ 

If  we  set  \* 
becomes, 


i  =  l,  2,  and  =  a£n[ 


(b+Xi )  (b+X'->)^ 

I  t - ^],  then  (2.5) 


b(b+X  j^-hX2) 


(2.6)  =  exp(-(XjTj  +  X^T2  +  ^ 

this  is  the  bivariate  exponential  distribution  (BVE)  of  Marshall  and  Olkin  (1967). 
The  marginal  distributions  corresponding  to  (2.6),  in  terms  of  the  parameters  a  and 

8,  are 


fb+X; 


(2.7) 


F-(t)  =  exp[  — a£n[— ^]Tj,  i  =  l,2,  and  r  ^  0. 
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If  we  let  o.i~)  =  o.T^,  for  constants  a,  S  >  0,  and  proceed  as  above,  then 
the  resulting  bivariate  survival  function  is  the  bivariate  ^eibull  of  Marshall  and 
Olkin  (196'/). 

Contrast  the  above  results  to  those  obtained  by  Lindley  and  Singpurwalla 
(1986),  who,  in  effect,  describe  their  uncertainty  about  Aj(t),  the  cumulative  rate  of 
component  i,  by  a  gamma  distribution,  and  obtain  the  logistic  and  the  Pareto  as  their 
bivariate  and  univariate  survival  functions,  respectively.  The  bivariate  logistic 
obtained  by  the  above  authors  is  absolutely  continuous,  whereas  the  bivariate 
exponential,  obtained  by  describing  A-Ct)  by  a  gamma  process,  has  a  singular 
compcnent. 

In  the  next  section,  we  explain  as  to  why  a  consideration  of  the  gamma  leads 
US  to  the  BV'E. 

3.  The  Gamma  Process  as  a  Poisson  Shock  Model  with  Varying  Shock  Intensiti^ 

The  BVE  has  been  motivated  [cf.  Barlow  and  Proschan  (1975),  p.  136]  via  a 
consideration  of  three  independent  Poisson  shock  processes,  with  the  shocks  within 
a  process  being  nonfatal  and  having  equal  intensities.  The  latter  assumption  implies 
that  the  probability  of  failure  due  to  each  shock  is  the  same  for  all  the  shocks 
within  a  process.  In  what  follows,  we  point  out  that  the  gamma  process  environment 
boils  down  to  a  consideration  of  a  nonfatal  Poisson  shock  process  with  the  shocks 
having  different  intensities;  that  is,  each  shock  induces  its  own  probability  of 
failure. 

To  ascertain  the  above,  it  is  useful  to  recall  that  a  process  with  independent 
increments  having  no  Gaussian  components  and  no  fixed  points  of  discontinuity 
such  as  a  gamma  process  -  can  be  represented  as  a  sum  of  a  countable  number  of 


jumps  of  random  height  at  a  countable  number  of  random  points  [cf.  Ferguson  and 


Klass  (1972)].  Furthermore,  if  for  any  fixed  e  >  0,  N(t,  e)  denotes  the  number  of 
jumps  in  (0,  t)  of  magnitude  ^  e,  then  {N(t,  e);  t  ^  0}  is  a  Poisson  process  with 


intensity  (M(e).  For  the  gamma  process,  M(e)  =  p 


u'^e*^  du,  where  p  is  a 


parameter  of  the  process;  see  for  example  Basawa  and  Prakasa  Rao  (1980),  p.  106. 

Let  U(e)  denote  the  magnitude  of  a  jump,  given  that  it  is  ^  e  >  0.  Then,  for  the 
gamma  process,  the  density  of  U(e)  is  a  truncated  (at  e)  gamma  density  with  shape  0. 


To  invoke  the  above  results  for  the  purpose  at  hand,  we  consider  jumps  of 
size  Ax,  and  remark  that  jumps  of  size  X,  where  X  ^  kAx,  for  k  =  l,  2,  .  .  .  ,  occur 
in  accordance  with  a  Poisson  process  with  intensity  M(kAx),  and  that  U(kAx)  is  a 
truncated  gamma  density  with  shape  0.  Thus  jumps  of  size  X,  where 
X  e  [kAx,  (k  +  l)Ax),  occur  in  accordance  with  a  Poisson  process  with  intensity 
AM  (kAx),  where  AM  (kAx)  =  M((k+l)Ax)  M(kAx),  k  =  l,  2,  .  .  .,  .  A  jump  in  the 

cumulative  failure  rate  can  be  viewed  as  the  result  of  a  discrete  shock,  the  effect 
of  which  IS  to  cause  the  failure  of  any  one,  or  both  components  of  the  system.  We 
would  expect  that  the  probability  of  failure  of  a  component  is  a  function  of  the 
size  of  the  shock,  and  if  Ax  is  sufficiently  small,  the  probability  of  failure  would 
be  approximately  the  same  for  all  X  G  [kAx,  (k+l)Ax). 


Let  PQj(kAx)  [pjQ(kAx)]  be  the  probability  that  a  shock  of  size 
X  £  [kAx,  (k+l)Ax)  will  cause  Component  1  [2]  to  fail  and  Component  2[1]  to 
survive.  Let  PqqILAx)  [p^j(kAx)]  be  the  probability  that  a  shock  will  cause  both 
components  to  fail  [survive].  Thus,  for  shocks  of  size  X  t  [kAx,  (k+l)Ax), 
k  =  l,  2,  .  .  .  ,  we  have  a  standard  nonfatal  shock  model  with  only  a  common  shock 
applied  to  both  components.  Thus,  for  any  t-,  ^  0,  the  probability  that  the 
components  will  survive  shocks  of  magnitude  X  £  [kAx,  (k+l)Ax)  is. 
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(3.1)  P(T^  >  Tp  T•^  >  T^:X  £  [kAx,  (k+l)Ax))  = 
exp{-AM(kAx)  pg^lkAxlr^  -  AM(kAx)  P^glkAxlr^  -  AM(kAx) 

max  (Tp  T-,)}. 

The  probability  that  the  system  will  survive  shocks  of  all  magnitudes  X  is 
the  product  of  (3T)  over  all  k.  Taking  the  limit  as  Ax  i  0,  we  have 

pCC 

(3.2)  P(T,  >  T.,  >  T.,)  =  exp{-  pQ^(x)dM(x)  - 

^  ^  -  Jq 

cc 

pQg(x)dM(x)  niax(Tp 

0 

and  this  is  of  the  form  (2.6). 

A  comparison  of  (3.2)  with  (2.6)  enables  us  to  interpret  Xp  and  X|^  in  term 
of  the  survival  probabilities  of  jumps  of  magnitude  X  and  their  arrival  intensities, 


p^Q(x)dM(x)T. 

0 


in  a  representation  of  the  gamma  process. 
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